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ABSTRACT 

Since the discovery of the first bona-fide brown dwarfs and extra-solar planets 
in 1995, the field of low mass stars and substellar objects has considerably 
progressed, both from theoretical and observational viewpoints. Recent 
developments in the physics entering the modeling of these objects have led 
to significant improvements in the theory and to a better understanding of 
their mechanical and thermal properties. This theory can now be confronted 
with observations directly in various observational diagrams (color-color, color- 
magnitude, mass-magnitude, mass-spectral type), a stringent and unavoidable 
constraint which became possible only recently, with the generation of synthetic 
spectra. In this paper, we present the current state-of-the-art general theory 
of low-mass stars and sub-stellar objects, from one solar mass to one Jupiter 
mass, regarding primarily their interior structure and evolution. This review is 
a natural complement to the previous review on the atmosphere of low-mass 
stars and brown dwarfs (Allard et al 1997). Special attention is devoted to the 
comparison of the theory with various available observations. The contribution 
of low-mass stellar and sub-stellar objects to the Galactic mass budget is also 
analysed. 
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1. INTRODUCTION 

Interest in the physics of objects at the bottom of and below the main sequence 
(MS) originated in the early demonstration by Kumar (1963) that hydrogen-burning in a 
stellar core no longer occurs below a certain mass and that below this limit hydrostatic 
equilibrium against gravitational collapse is provided by electron degeneracy. Simple 
analytical arguments, based on the balance between the classical ionic thermal pressure 
and the quantum electronic pressure yield for this H-burning minimum mass (HBMM) 
tohbmm ~ 0.1 Mq, while the first detailed evolutionary calculations gave ttihbmm ~ 0.085 M 
(Grossman et al 1974). Tarter (1975) proposed the term "brown dwarfs" (BD) for objects 
below this H-burning limit. D'Antona and Mazzitelli (1985) first pointed out that the 
luminosity below wihbmm would stretch by a few orders of magnitude over a few hundredths 
of a solar mass, making the observation of BDs a tremendously difficult task. Subsequent 
benchmarks in low-mass star (LMS) and BD theory are due primarily to VandenBerg et al 
(1983), D'Antona and Mazzitelli (1985), Dorman et al (1989), Nelson et al (1986, 1993a) 
and the Tucson group (Lunine et al 1986, 1989, Burrows et al 1989, 1993). In spite of this 
substantial theoretical progress, all of these models failed to reproduce the observations 
at the bottom of the MS and thus could not provide a reliable determination of the 
characteristic properties (mass, age, effective temperature, luminosity) of low-mass stellar 
and substellar objects. A noticeable breakthrough came from the first calculations by a few 
groups of synthetic spectra and atmosphere models characteristic of cool (T cS < 4000 K) 
objects (Allard 1990, Brett & Plez 1993, Saumon et al 1994, Allard & Hauschildt 1995, Brett 
1995, Tsuji et al 1996, Hauschildt et al 1999). This allowed the computation of consistent 
non-grey evolutionary models (Saumon et al 1994, Baraffe et al 1995) and the direct 
confrontation of theory and observation in photometric passbands and color- magnitude 
diagrams, thus avoiding dubious transformations of observations into theoretical L-T e ff 
Hertzsprung-Russell diagrams. 
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In the meantime, the search for faint (sub)stellar objects has bloomed over the past 
few years. Several BDs have now been identified since the first discoveries of bona-fide BDs 
(Rebolo et al 1995, Oppenheimer et al 1995) either in young clusters (see Basri, this volume 
and Martin 1999 for reviews, and references therein) or in the Galactic field (Ruiz et al 1997, 
Delfosse et al 1997, Kirkpatrick et al 1999a). The steadily increasing number of identified 
extra-solar giant planets (EGPs) since the discovery of 51PegB (Mayor & Queloz 1995, 
see Marcy & Butler 1998 for a review) has opened up a new era in astronomy. Ongoing 
and future ground-based and space-based optical and infrared surveys of unprecedented 
faintness and precision are likely to reveal hundreds more of red dwarfs, brown dwarfs and 
giant planets which will necessitate the best possible theoretical foundation. The correct 
understanding of the physical properties of these objects bears major consequences for a 
wide range of domains of physics and astrophysics: dense matter physics, planet and star 
formation and evolution, galactic evolution, missing mass. 

A general outline of the basic physics entering the structure and the evolution of BDs 
can be found in the excellent reviews of Stevenson (1991) and Burrows & Liebert (1993). 
It is the aim of this review to summarize the most recent progress realized in the theory of 
LMS, BDs and EGPs. We will focus mainly on the internal structure and the evolution 
of these objects since a comprehensive review on the atmosphere of LMS and BDs has 
appeared recently (Allard et al 1997). We will also consider the implications of LMS and 
BDs in a more general galactic context and evaluate their contribution to the Galactic mass 
budget. 



2. THE PHYSICS OF DENSE OBJECTS 



2.1. Interior Physics 
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2.1.1. Equation of State 

Central conditions for LMS, hereafter identified generically as objects below a solar 
mass, and substellar objects (SSO) for solar composition range from a maximum density 
p c ~ 10 3 gcm~ 3 at the hydrogen-burning limit (m 0.07 M , see below) to p c ~ 10 gem -3 
for Saturn (= 5 x 1O~ 4 M ) at 5 Gyr, and from T c ~ 10 7 K for the Sun to T c ~ 10 4 
K for Saturn at the same age, spanning several orders of magnitudes in mass, density 
and temperature. Effective temperatures range from ~ 6000 K to ~ 2000 K in the 
stellar domain and extend down to ~ 100 K for Saturn. Molecular hydrogen and other 
molecules become stable for kT <3 x 10~ 2 Ryd (Tsi5x 10 3 K), a condition encountered 
in the atmosphere, or even in deeper layers, of most of these objects. Under these 
conditions, the interior of LMS and SSO is essentially a fully ionized H + /He ++ plasma 
characterized by coupling parameters Tj = (Z 5 ^ 3 )e 2 /a e kT rs 0.1-50 for classical ions and 
r s = (Zy^a/ao « 0.1-1 for degenerate electrons. With I\ = (Z 5 / 3 )(2.693 x 10 5 K/T)n 2 { 3 , 
r s = 1.39/(p//i e ) 1 / 3 = 1.172 n,2^ 3 , n 2 A = n e /10 24 cm -3 ^ (p/1.6605 gem -3 )//" 1 the electron 
number-density, a = (£:^-) 1//3 = a e(-^) 1 ^ 3 the mean interionic distance, /i -1 = (Z) / (A) the 
electron mean molecular weight, p the mass-density, A the atomic mass, a the electronic 
Bohr radius and <> denotes the average number- fraction (X) = J2iXiXi (x, = Ni/ ^Ni). 
The temperature is of the order of the electron Fermi temperature Tp so that the degeneracy 
parameter ip = kT/kT F w 3.314 x 10 -6 T (p e / p) 2/3 is of the order of unity. The classical 
(Maxwell-Boltzman) limit corresponds to ip — > +oo, whereas ip — > corresponds to 
complete degeneracy. The afore-mentioned thermodynamic conditions yield ip ~ 2 — 0.05 in 
the interior of LMS and BDs implying that finite-temperature effects for the electrons must 
be included to describe accurately the thermodynamic properties of the correlated, partially 
degenerate electron gas. Moreover, the Thomas-Fermi wavelength \ T f = {kT F / {Snriee 2 )) 1 ! 2 
is of the order of a, so that the electron gas is polarized by the external ionic field and 
electron-ion coupling must be taken into account in the plasma hamiltonian. Last but 
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not least, the electron average binding energy can be of the order of the Fermi energy 
Ze 2 /ao ~ ep so that presswre-ionization takes place along the internal profile. Figure [3] 
illustrates central characteristic quantities for LMS and SSOs from the Sun to Jupiter. 

Above ~ 0.4 M & , the structure evolves slowly with increasing mass from a n=3/2 
towards a n = 3 polytrope, which yields the correct P c for the Sun, due to the growing 
central radiative core. This leads to increasing central pressures and densities for increasing 
mass (increasing polytropic index) in this mass range. Below ~ 0.3-0.4 M , the core 
becomes entirely convective and follows the behaviour of a n = 3/2 polytrope. Since the 
gas is still in the classical regime (ip > 1), m oc R and the central density increases with 
decreasing mass, p c oc m~ 2 . Below about the H-burning limit, electron degeneracy becomes 
dominant (ip<0.1), so that one approaches the relation m oc R~ 3 (for ip=0) and density 
decreases again with decreasing mass, p c oc m 2 . These various effects yield a non-monotomic 
behaviour of the central density and pressure with mass with a minimum in the stellar 
regime around 0.4 M & and a maximum near the H-burning limit. 

The equation of state (EOS) of low-mass objects thus requires a detailed description 
of strongly correlated, polarisable, partially degenerate classical and quantum plasmas, 
plus an accurate treatment of pressure partial ionization, a severe challenge for dense 
matter physicists. Several steps towards the derivation of such an accurate EOS for dense 
astrophysical objects have been realized since the pioneering work of Salpeter (1961), with 
major contributions by Fontaine et al (1977), Magni and Mazzitelli (1979), Marley and 
Hubbard (for objects with m < 0.2 M Q ), Hummer & Mihalas (1988) (primarily devoted to 
conditions characteristic of solar-type stellar envelopes) and Saumon and Chabrier (Chabrier 
1990, Saumon & Chabrier 1991, 1992, Saumon et al 1995, hereafter SCVH). We refer the 
reader to Saumon (1994) and SCVH for a detailed review and an extensive comparison 
of these different EOS in the domain of interest. The SCVH EOS presents a consistent 
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treatment of pressure ionization and includes significant improvements with respect to 
previous calculations in the treatment of the correlations in dense plasmas. It compares 
well with available high-pressure shock-wave experiments in the molecular domain, and 
with Monte-Carlo simulations in the fully-ionized, metallic domain (see references above 
for details). Recently laser-driven shock- wave experiments on D 2 have been conducted 
at Livermore (Da Silva et al 1997, Collins et al 1998) which reach for the first time the 
pressure-dissociation and ionization domain and thus probe directly the thermodynamic 
properties of dense hydrogen under conditions characteristic of BDs and giant planets (GP). 
The relevance of these experiments for the interior of SSOs can be grasped from Figure 
1 of Collins et al (1998). These experiments have shown the good agreement between 
the predictions of the Saumon-Chabrier model and the data (although there is certainly 
room for improvement). In particular the strong compression factor arising from hydrogen 
pressure-dissociation and ionization observed in the experiment (p/pi ~ 5.8) agrees well 
with the predicted theoretical value (see Figure 3 of Collins et al 1998). Any EOS devoted 
to the description of the interior of dense astrophysical objects must now be confronted 
with these available data. 

The SCVH EOS is a pure hydrogen and helium EOS, based on the so-called 
additive- volume-law between the pure components (H and He). The accuracy of the 
additive- volume-law has been examined in detail by Fontaine et al (1977). The invalidity 
of this approximation to describe accurately the thermodynamic properties of the mixture 
is significant only in the partial ionization region (see SCVH), which concerns only a 
few percent of the stellar mass under LMS and BD conditions. The effect of metals 
on the structure and the evolution of these objects has been examined in detail by 
Chabrier & Baraffe (1997, §2.1). As shown by these authors, because of their negligible 
number-abundance (~ 0.2%), metals do not contribute significantly to the EOS and barely 
modify the structure and evolution (~ 1% in T eff and ~ 4% in L) of these objects. In the 
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low-density limit characteristic of the atmosphere, the SCVH EOS recovers the perfect gas 
limit and thermal contributions from various atomic or molecular species can be added 
within the afore-mentioned additive-volume law formalism, which is exact in this regime. 
Only in the (denser) metal-depleted atmospheres (Z < 1CT 2 Z Q ) of the densest objects 
(<O.1M ), one finds a slight departure from ideality, with a 1% to 4% effect on the 
adiabatic gradient V ' a d (Chabrier & Baraffe 1997). 

2.1.2. Nuclear rates. Screening factors 

Although the complete PP chain is important for nucleosynthesis, the thermonuclear 
processes relevant from the energetic viewpoint under LMS and BD conditions are given by 
the PPI chain (see e.g. Burrows & Liebert 1993, Chabrier & Baraffe 1997) : 

p + p -> al + e + + u e ; p + d -> 3 He + 7; 3 He + 3 He ^ 4 He + 2p (1) 

Below ~ O.7M , the PPI chain contributes to more than 99% of the energy generation 
on the zero-age Main Sequence (ZAMS), and the PPII chain to less than 1%. The 
destruction of 3 He by the above reaction is important only for T > 6 x 10 6 K, i.e masses 
m>O.15M for ages t < 10 Gyr, for which the lifetime of this isotope against destruction 
becomes eventually smaller than a few Gyr. 

As examined in §4.6, the abundance of light elements (D, Li, Be, B) provides a powerful 
diagnostic to identify the age and/or the mass of SSOs. The rates for these reactions (see 
e.g. Nelson et al 1993b) in the vacuum, or in an almost perfect gas where kinetic energy 
largely dominates the interaction energy, are given by Caughlan and Fowler (1988) and 
Ushomirsky et al (1998) for updated values. The reaction rate R (in cm _3 s _1 ) in the 
vacuum is given by Rq oc e~ 3e0//fcT where eo corresponds to the Gamow-peak energy for 



- 12 - 



non- resonant reactions (Clayton 1968). However, as mentioned in the previous section, 
non-ideal effects dominate in the interior of LMS and BDs and lead to polarization of the 
ionic and electronic fluids. These polarization effects due to the surrounding particles yield 
an enhancement of the reaction rate, as first recognized by Schatzman (1948) and Salpeter 
(1954). The distribution of particles in the plasma reads : 



n(r) = ne- Ze<t>ir)/kT , with Mr) = — + ^(r) (2) 

r 

where ip(r) is the induced mean field potential due to the polarization of the surrounding 
particles. This induced potential lowers the Coulomb barrier between the fusing particles 
and thus yields an enhanced rate in the plasma R— E x R where 



E = lim r ^ Q {g 12 (r)exp( 1 *f )} (3) 

is the enhancement (screening) factor and g 12(f) the pair-distribution function of particles 
in the plasma. Under BD conditions, not only ion screening must be considered but also 
electron screening, i.e. E = E,- t x E e . Both effects are of the same order (E$ ~ E e ~ a few) 
and must be included in the calculations for a correct estimate of the light element-depletion 
factor (see §4.6). Figure |2| portrays the evolution of the central temperature for objects 
respectively above, at the limit of and below the hydrogen-burning minimum mass, with 
lines indicating the hydrogen, lithium and deuterium burning temperatures in the plasma. 



2.1.3. Transport Properties 



Energy in the interior of LMS below ~ 0.4 M , BDs and GPs is transported essentially 
by convection (see e.g. Stevenson 1991). According to the mixing length theory (MLT), the 
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convective flux reads (Cox and Giuli 1968) : 

Fconv OC P V conv C P 5T OC (gV2 /r V2 )( J_ )2pCpCsT(v _ ^3/2 (4) 

Hp 

where Q = — (ffffOp is the volume expansion coefficient, Cp is the specific heat at constant 
pressure, cs = (TiP/ p) 1 / 2 is the adiabatic speed of sound, 5T is the temperature excess 
between the convective eddy and the surrounding ambient medium, V the temperature 
gradient, f con v is the convective velocity, /andifp denote the mixing length and the pressure 
scale height, respectively, and V e ~ V a d is the eddy temperature gradient. The last term 
in parenthesis on the r.h.s. of (4) defines the fluid superadiabaticity, i.e. the fractional 
amount by which the real temperature gradient exceeds the adiabatic temperature gradient. 
Below a certain mass, the inner radiative core vanishes and the star becomes entirely 
convective (VandenBerg et al 1983, D'Antona & Mazzitelli 1985, Dorman et al 1989). This 
minimum mass, calculated with consistent non-grey atmosphere models and with the most 
recent OPAL radiative opacities (Iglesias and Rogers 1996) for the interior, is found to be 
n^conv — 0.35 M within the metallicity range 10~ 2 < Z/Z Q < 1 (Chabrier & Baraffe 1997 
§3.2). The Rayleigh number is defined as Ra = gQ ^ p ( V — Vad), where g ~ 10 3 -10 5 cm s~ 2 
is the surface gravity, £ is the thermal diffusivity (conductive or radiative diffusivity) and v 
is the kinematic viscosity. In the interior of LMS, Ra ~ 10 25 so that convection is almost 
perfectly adiabatic and the MLT provides a fairly reasonable description of this transport 
mechanism. Variation of the MLT parameter a = I /Hp between 1 and 2 in the interior is 
found to be inconsequential below ~ 0.6 M Q (see Baraffe et al 1997). 

However convection can be affected or even inhibited by various mechanisms. Maximum 
rotation velocities for LMS and BDs are usually in the range ~ 20 — SOkms" 1 (see e.g 
Delfosse et al 1998a) with values as large as 50 kms" 1 and 80 kms -1 in the extreme case 
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of Kelu 1 (Ruiz et al 1997, Basri et al 2000). This corresponds to an angular velocity 
0~5x 10~ 4 rad s _1 for a characteristic radius R ~ 0.1 R e (£lj up = 1.76 x 10 -4 rad s _1 ). 
Within most (> 95% in mass) of the interior of LMS and BDs, t> CO nv is of the order of 
10 2 cm s^ 1 (<< cs) and H P ~ 10 9 cm. Thus, the Rossby number Ro = v conv /ill rs 10~ 3 , 
where / ~ H P , and convection can in principle be inhibited by rotation. However, the fact 
that lithium is not observed in objects with m>0.06M Q (see §4.6) suggests that some 
macroscopic transport mechanism (e.g. meridional circulation, turbulence or convection) 
remains efficient throughout the star. The Reynolds number Re = v conv l/u remains largely 
above unity in stellar and SSO interiors so that convection is not inhibited by viscosity. 
Magnetic inhibition requires a magnetic velocity va = (B 2 / Anp) 1 / 2 > v conv (Stevenson 1979, 
1991), i.e. B > 10 4 Gauss under the conditions of interest, only slightly above the predicted 
fields in these objects (see §2.3). However, it has been proposed by Stevenson (1979) 
that in a rapidly rotating fluid with a significant magnetic field, the Proudman- Taylor 
theorem no longer applies so that vertical convective motions could be possible. Lastly, 
convection efficiency can be diminished by the presence of a density gradient. As shown by 
Guillot (1995), this may occur in the interior of gaseous planets, due to a gradient of heavy 
elements, leading to inefficient convective transport in the envelope of these objects. 

The treatment of convection in the outer (molecular) layers, above and near the 
photosphere, is a more delicate question. The Rayleigh number in this region is Ra ~ 10 15 , 
a rather modest value by the usual standards in turbulence. Since the MLT, by definition, 
applies to the asymptotic regime Ra — > oo, it is no longer valid near the photosphere. 
Much work has been devoted to the improvement of this formalism and to the derivation 
of a non-local treatment of convection. One of the most detailed attempts is due to 
Canuto and Mazitelli (Canuto & Mazzitelli 1991, CM). The original CM formalism was 
based on a linear stability analysis, whereas energy transport by turbulence is a strongly 
non-linear process. It has been improved recently by including to some extent non-linear 
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modes in the energy rate (Canuto et al 1996) but it yields similar results as the initial 
CM model. The formalism, however, still requires the calibration of a free parameter, 
which represents a characteristic "mixing" scale and in this sense resembles the MLT. A 
more severe shortcoming of the CM formalism is that the predicted outermost limit of 
the convection zone for the Sun and the amount of superadiabaticity in this zone are in 
significant disagreement, both quantitatively and qualitatively, with the results obtained 
from coupled hydrodynamic-radiation 3D simulations (Demarque et al 1999, Nordlund & 
Stein 1999). The 3D simulations yield excellent agreement with the observational data at 
several diagnostic levels: the thermal structure (and the inferred depth of the convection 
zone), the dynamical structure (vertical velocity amplitude and spectral line synthesis) and 
the p-mode frequencies and amplitudes for the Sun, without the use of free parameters (see 
e.g. Spruit & Nordlund 1990 for a review, and Nordlund & Stein 1999). This provides now 
a high degree of confidence in 3-D hydrodynamic models of stellar surface layers and the 
inferred transition from convective to radiative energy transport. 

The thermal structure, for example, is a very robust property of the numerical models, 
since it depends relatively little on the level of turbulence (and thus on numerical resolution) 
(Stein & Nordlund 2000; Nordlund & Stein 1999). These results illustrate the fact that 
any formalism based on a homogeneous description of convection (e.g. MLT, CM) can 
not describe accurately this strongly inhomogeneous process. Interestingly enough, the 
standard MLT is found to compare reasonably well with these simulations, at least for 
the thermal profile (see afore- mentioned references), and thus seems to offer a reasonable 
(or least worst!) overall description of convective transport, even in the small-efficiency 
convective regions. Indeed, in contrast to 1-D models computed with the CM formalism, 
the 3-D models do not have a steeper and more narrow superadiabatic structure for the 
Sun than obtained with the classical Bohm-Vitense mixing length recipe (Nordlund & 
Stein 1999). It it thus fair to say that, in the absence of a correct non-local treatment of 
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convection in LMS and BD interiors, the standard MLT is probably the most reasonable 
choice, at least for the present objects, LMS and SSOs. Clearly the developpement of 3-D 
hydrodynamics models of the atmosphere of these objects, and the calibration of the mixing 
length from these simulations, as done for solar-type stars (Ludwig et al 1999), represents 
one of the next major challenges in the theory. 

Another possible mechanism of energy transport in stellar interiors is conduction. Its 
efficiency can be estimated as follows: the distance / over which the temperature changes 
significantly is I ~ (x^) 1//2 > where \ * s the thermal diffusivity and t is the time during 
which the temperature change occurs. For x ~ 10" 1 cm 2 s _1 , characteristic of metals, and 
t ~ 10 9 yr, I ~ 10 2 km. Heat can thus possibly be transported by conduction only over 
a limited range of the interior, providing the density is high enough and the temperature 
low enough for electron conductivity to become important. Indeed, below the HBMM 
~ 0.07 Mq, the interior becomes degenerate enough during cooling that the conductive 
flux F C ond = -KcondVT, where A" CO nd = T3 and /t con d is the conductive opacity, 
becomes larger than the convective flux. Old enough BDs in the mass-range 0.02-0.07 M & 
become degenerate enough to develop a conductive core, which slows down the cooling L(t) 
(Chabrier et al 2000). 



2.2. Atmosphere 

Knowledge of the atmosphere is needed for two reasons : (i) as a boundary condition 
for the interior profile in the optically-thick region and (ii) as a description of the emergent 
radiative flux. A comprehensive review of the physics of the atmosphere of LMS and BDs 
can be found in Allard et al (1997). Only a general outline of the main characteristics of 
these atmospheres will be mentioned in the present review. 



2.2.1. Spectral distribution 



Surface gravity g = Gm/R 2 for LMS and SSOs range from log g ~ 4.4 for a 
main sequence 1.0 M star to logg ~ 3.4 for Jupiter, with a maximum logg ~ 5.5 at 
the H-burning limit, for solar metallicity. This yields P p h ~ g/R ~ 0.1 — 10 bar and 
p p h ~ 10~ 6 -10~ 4 gcur 3 at the photosphere. Collision effects become significant under 
these conditions and induce molecular dipoles on e.g. H 2 or He-H 2 , yielding so-called 
collision-induced absorption (CIA) between roto-vibrational states of molecules which 
otherwise would have only quadrupolar transitions (Linsky 1969; Borysow et al 1997). At 
first order this CIA coefficient scales as kcia oc n 2 H2 x kh 2 h 2 an d thus becomes increasingly 
important as soon as H 2 molecules become stable. The CIA of H 2 suppresses the flux 
longward of 2 /zm in the atmosphere of LMS, BDs and GPs. This is one of the reasons (the 
main one for metal-depleted objects) for the redistribution of the emergent radiative flux 
toward shorter wavelengths in these objects (Lenzuni, Chernoff & Salpeter 1991; Saumon 
et al 1994, Allard & Hauschildt 1995, Baraffe et al 1997). The CIA of H 2 and H 2 -He, and 
the bound-free and free-free opacities of H and provide the main continuum opacity 
sources below ~ 5000 K for objects with metal-depleted abundances. Below T eff < 4000 
K, most of the hydrogen is locked into H 2 and most of the carbon in CO. Excess oxygen 
is bound in molecules such as TiO, VO and H 2 0, with some amounts also in OH and O 
(see e.g. Figure 1 of Fegley & Lodders 1996). Metal oxides and metal hydrides (FeH, 
CaH, MgH) are also present. The energy distribution of solar-abundance M-dwarfs is thus 
entirely governed by the line absorption of TiO and VO in the optical, and H 2 and CO 
in the infrared, with no space for a true continuum. The VO and TiO band strength index 
are used to classify M-dwarf spectral types (see e.g. Kirkpatrick et al 1991). The strongly 
frequency-dependent absorption coefficient due to the line transitions of these molecules, as 
well as Rayleigh scattering (oc z/ 4 ) shortward of ~ 0.4 /im, yield a strong departure from a 
black-body energy distribution (see e.g. Figure 5 of Allard et al 1997). 
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At T cS 2000 K, near the H-burning limit, signatures of metal oxides and hydrides 
(TiO, VO, FeH, CaH bands) disappear from the spectral distribution (although some 
amount of TiO remains present in the atmosphere), as observed e.g. in GD 165B 
(Kirkpatrick et al 1999b). Alkalies are present under their atomic form. The disappearance 
of TiO-bands prompted astronomers to suggest a new spectral type classification, the 
"L"-type, for objects with T c g below the afore-mentioned limit (Martin et al 1997, 1998, 
Kirkpatrick et al 1999a). Below a local temperature T ~ 1300 - 1500 K for P ~ 3 - 10 
bars, for a solar abundance-distribution, carbon monoxide CO was predicted to dissociate 
and the dominant equilibrium form of carbon was predicted to become CH 4 (Allard & 
Hauschildt 1995, Tsuji et al 1995, Fegley & Lodders 1996). Note that this transition occurs 
gradually with some of the two elements being present in the stability field of the other (see 
e.g. Fegley & Lodders 1996, Lodders 1999) so that cool stars could contain some limited 
amount of methane and CO may be visible in objects with T e g < 1800 K, as indeed detected 
in G1229B (Noll et al 1997; Oppenheimer et al 1998). This prediction has been confirmed 
by the spectroscopic observation of Gliese229B and the identification of methane absorption 
features at 1.7, 2.4 and 3.3 /iin (Oppenheimer et al 1995). The presence of methane in its 
spectrum confirmed unambiguously its sub-stellar nature (Oppenheimer et al 1995, Allard 
et al 1996, Marley et al 1996). Objects like G1229B, characterized by the strong signature of 
methane absorption in their spectra, define a new spectral class of objects, the "methane" 
(sometimes called "T") BDs. Methane absorption in the H and K bands yields a steep 
spectrum at shorter wavelengths and thus blue near-infrared colors, with J — K < but 
I — J>5 (Kirkpatrick et al 1999a, Allard 1999). Although the transition temperature 
between "L" and "methane" dwarfs is not determined precisely yet, it should lie in the 
range 1000<T cfr < 1700 K (see §4.5.2). 

Below T c g ~ 2800 K, complex O-rich compounds condense in the atmosphere, slightly 
increasing the carbomoxygen abundance ratio (see e.g. Tsuji et al 1996, Fegley & Lodders 
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1996, Allard et al 1997). Different constituents will condense at a certain location in the 
atmosphere with an abundance determined by chemical equilibrium conditions (although 
non-equilibrium material may form, depending on the time scale of the reactions, as 
mentioned below) between the gas phase and the condensed species. The formation of 
condensed species depletes the gas phase of a number of molecular species (e.g. VO, TiO 
which will be sequestered into perovskite CaTiOa), modifying significantly the emergent 
spectrum (see e.g. Fegley & Lodders 1994 for the condensation chemistry of refractory 
elements in Jupiter and Saturn). The equilibrium abundances can be determined from the 
Gibbs energies of formation either by minimization of the total Gibbs energy of the system 
(Sharp & Huebner 1990, Burrows & Sharp 1999), or by computing the equilibrium pressures 
of each grain species (Grossman 1972, Alexander & Ferguson 1994, Allard et al 1998). At 
each temperature, the fictitious pressure Pc of each condensed phase under consideration 
is calculated from the partial pressures of the species % which form the condensate (e.g. 
Al and O for corundum Al 2 0^) Pi = N{IZT , where iVj is the number of moles of species 
i and 1Z is the gas constant, determined by the vapor phase equilibria. This fictitious 
pressure Pc is compared to the equilibrium pressure P eq , calculated from the Gibbs energy 
of formation of the condensate. The abundance of a condensed species is determined by 
the condition that this species be in equilibrium with the surrounding gas phase, Pc > P eq 
(Grossman 1972). The opacities of the grains are calculated from Mie theory. As suggested 
by Fegley & Lodders (1996) and Allard et al (1998), refractory elements, Al, Ca, Ti, Fe and 
V are removed from the gaseous atmosphere by grain condensation at about the corundum 
(AI2O3) or perovskite (CaTi0 3 ) condensation temperature T< 1800 K. Rock-forming 
elements (Mg, Si, Fe) condense as iron and forsterite (Mg 2 Si0 4 ) or enstatite (MgSi0 3 ) 
within about the same temperature range (depending on P). Therefore the spectral features 
of all these elements will disappear gradually for objects with T e s below these temperatures 
(see e.g. Fegley & Lodders 1994, 1996, Lodders 1999, Burrows & Sharp 1999 for detailed 
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calculations). For jovian- type effective temperatures (T eS ~ 125 K), H 2 and NH 3 condense 
near and below the photosphere, and water and ammonia bands disappear completely for 
T eff < 150 K and 80 K, respectively (Guillot 1999). 

The gas abundance strongly depends on pressure and temperature so that the 
abundances of various species varies significantly with the mass (and T eS ) of the 
astrophysical body. As shown e.g. in figure 2 of Lodders (1999), the condensation point of 
the dominant clouds lies much closer to the photosphere for M-dwarfs than for G1229B or 
- worse - for Jupiter. In other words the location of a given grain condensation lies deeper 
in Jupiter than in G1229B. For late M-dwarfs and for massive and/or young BDs, the main 
cloud formation is predicted to occur very near the photosphere. This is consistent with the 
fact that all the DENIS and 2MASS objects discovered near the bottom of and below the 
MS exhibit strong thermal heating and very red colors (Delfosse et al 1998b, Kirkpatrick et 
al 1999a) (see §4.5.2). Indeed the atmospheric heating due to the large grain opacity (the 
so-called greenhouse or backwarming effect), and the resulting enhanced H 2 dissociation, 
yield a redistribution of the IR flux, as proposed initially by Tsuji et al (1996). 

A key question for the grain formation process is the size of the grains (see e.g. 
Alexander & Fergusson 1994). The suppression of the flux in the optical in G1229B suggests 
grain sizes of ~ 0.1 /im as a source of continuum opacity (Jones & Tsuji 1997, Griffith et al 
1998), although the wings of alkali resonance lines (K I, Na I) is also a source of absorption 
in this region (Tsuji et al 1999, Burrows et al 2000). Recent calculations assume a grain-size 
distribution in the submicron range (Allard et al 1999). Inclusion of the grain absorption 
with this size distribution in the atmosphere of objects near the limit of the MS successfully 
reproduces the observed colors of GD 165B, Kelu-1 and of DENIS objects (Leggett et al 
1998, Allard 1999, Goldman et al 1999, Chabrier et al 2000). However, the spectrum of 
GL229B shows no indication for dust in its IR spectrum from 1 to 5 jim (Allard et al 1996, 
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Marley et al 1996, Tsuji et al 1996, Oppenheimer et al 1998, Schultz et al 1998). This 
suggests ongoing dynamical processes such as grain settling in SSO atmospheres. Indeed, as 
noted by Chabrier et al (2000), although convection occurs only in optically-thick regions 
(r > 1) for L-dwarfs, the top of the convection zone lies only about one pressure scale height 
or less from the photosphere (r ~ 1) near which most grains condense for objects near the 
bottom of the MS. This can have important consequences on the formation and settling of 
atmospheric grains. Indeed, although the temperature at the top of the convective zone is 
found to be generally above the condensation temperature of all grains, convection-induced 
advection or turbulent diffusion could efficiently bring material upward to the region of 
condensation and maintain small-grain layers, which otherwise will settle gravitationally, 

In general, the grain formation process involves a balance between various dynamical 
timescales, such as the condensation, evaporation, coagulation, coalescence and convection 
timescales (see e.g. Rossow 1978), not mentioning large-scale hydrodynamical instabilities, 
typical of weather conditions on Earth. See for example Marley et al (1999) for an early 
attempt at modelling BD and EGP atmospheres. Moreover, although to an excellent 
first approximation the atmosphere of the solar jovian planets are in thermodynamical 
equilibrium, non-equilibrium species with chemical equilibrium timescales larger than the 
convection timescale can be dredged up to the photosphere by convection, as e.g. CO, 
PH 3 , GeH 4 in Jupiter (Fegley & Lodders 1994). A correct understanding and a reliable 
description of this complex grain formation process represent a major challenge for theorists 
and the observation of spectra of SSOs over a significant temperature-range is necessary to 
provide guidance through this maze. 
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2.2.2. Transport Properties 

As mentioned previously, below ~ 5000 K, i.e. m ~ 0.6-0.8 M Q depending on the 
metallicity, hydrogen atoms recombine, n# 2 increases, as does the opacity through H 2 
CIA. The radiative opacity k increases by several orders of magnitude over a factor 2 in 
temperature, which in turn decreases the radiative transport efficiency (Trad oc ^). On the 
other hand, the presence of molecules increases the number of internal degrees of freedom 
(vibration, rotation) and thus the molar specific heat C p , which in turn decreases the 
adiabatic gradient (dT/dP) ad = -^-p- Both effects strongly favor the onset of convection in 
the optically-thin (r < 1) atmospheric layers (Copeland et al 1970, Allard 1990, Saumon 
et al 1994, Baraffe et al 1995) with a maximum radial extension for the convection zone 
around T ef j ~ 3000 K. The combination of relatively large densities, opacities and specific 
heat in layers where H 2 molecules are present contribute to a very efficient convective 
transport. Convection is found to be adiabatic almost up to the top of the convective region 
and the extension of the superadiabatic layers is very small compared to solar type stars 
(Brett 1995, Allard et al 1997). Thus a variation of the mixing length between H P and 2H P 
barely affects the thermal atmospheric profile (Baraffe et al 1997 §3). 

The presence of convection in the optically thin layers precludes radiative equilibrium 
in the atmosphere (Tot ^ Tad) and requires the resolution of the transfer equation 
for radiative-convective atmosphere models for objects below ~ 0.7 M Q . This, and the 
afore-mentioned strong frequency-dependence of the molecular absorption coefficients, 
exclude grey model atmosphere and/or the use of radiative T(r) relationships to determine 
the outer boundary condition to the interior structure. As demonstrated by Saumon et al 
(1994) and Baraffe et al (1995, 1997, 1998), such outer boundary conditions overestimate 
the effective temperature and the HBMM and yield erroneous m-T eS and m-L relationships, 
two key relations for the calibration of the temperature-scale and for the derivation of 
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mass functions (see Chabrier & Baraffe 1997, §2.5 for a complete discussion). Correct 
evolutionary models for low mass objects require (i) the connection between the non-grey 
atmospheric (P,T) profile, characterized by a given (log g,T e s) and the interior (P,T) 
profile at a given optical depth r, preferably large enough for the atmospheric profile to be 
adiabatic and (ii) consistency between the atmospheric profiles and the synthetic spectra 
used to determined magnitudes and colors. For a fixed mass and composition, only one 
atmosphere profile matches the interior profile for the afore-mentioned boundary condition. 
This determines the complete stellar model (mass, radius, luminosity, effective temperature, 
colors) for this mass and composition (Chabrier & Baraffe 1997, Burrows et al 1997 for 
objects below 1300 K (their figure 1)). 

2.3. Activity 

The Einstein and Rosat surveys of the solar neighborhood and of several open clusters 
have allowed the identification of numerous M-dwarfs as faint X-ray sources, with typical 
luminosities of the order of the average solar luminosity L x ~ 10 27 erg s" 1 , reaching up to 
~ 10 29 erg s _1 in young clusters (see e.g. Randich, 1999 and references therein for a recent 
review). This suggests that late M-dwarfs are as efficient coronal emitters as other cool 
stars in terms of L x /L bo i, which expresses the level of X-ray activity. Moreover about 60% 
of M-dwarfs with spectral-type > M5 (Mboi — 12) show significant chromospheric activity 
with log(LHa/Lboi) — —3.9 (Hawley et al 1996). These observations provide important 
empirical relationships between age, rotation and activity. The L x /L bo i ratio seems to 
saturate above a rotational velocity threshold at a limit \og(L x / L bol ) ~ —3, implying that 
the intrinsic coronal emission of LMS and BDs is quite low, decreases with Ly and thus 
mass and does not increase with increasing rotation above the threshold limit. This suggests 
a saturation relation in terms of a Rossby number, i.e. the ratio of the rotational period over 
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the convective turnover time Ro = j^— (see §2.1.3). Because t conv increases for lower-mass 

tconv 

stars, they will saturate at progressively larger periods, i.e. lower rotational velocities. The 
velocity threshold for a 0.4 M star is estimated around 5-6 km s _1 (Stauffer et al 1997). 
Comparison of the Pleiades and the Hyades samples shows a steep decay of the X-ray 
activity in the Hyades. This result is at odds with a simple, monotonic rotation-activity 
relationship. The spin-down timescales for M-dwarfs are much longer than for solar-type 
stars. Hyades M-dwarfs do indeed show moderate or even rapid rotation (vsmi > 10 km 
s— 1) and thus should show strong X-ray emission. The fainter L x in the Hyades than in 
the Pleiades thus reinforces the saturation scenario. 

Delfosse et al (1998a) obtained projected rotational velocities and fluxes in the H a 
and H/3 lines for a volume-limited sample of 118 field M-dwarfs with spectral-type M0-M6. 
They found a strong correlation between rotation and both spectral type (measured by 
R — I colors) and kinematic population: all stars with measurable rotation are later than 
M3.5 and have kinematic properties typical of the young disk population, or intermediate 
between young disk and old disk. They interpret this correlation as evidence for a spin-down 
timescale that increases with decreasing mass, with this timescale being a significant 
fraction of the age of the young disk (~ 3 Gyr) at spectral type M4 (~ 0.15 M ). These 
data confirm the saturation relation inferred previously for younger or more massive stars: 
Lx/Lboi and L Ha /L bol both correlate with v sini for v sin i< 4-5 km s _1 and then saturate 
at 10~ 2 ' 5 and 10~ 35 , respectively. Recently, Neuhauser and Comeron (1998) have detected 
X-ray activity in young BDs with m ~ 0.04 M® in the Chamaleon star-forming region, with 
\ogLx = 28.41 and log Lx/L^i = —3.44. Coronal activity in BDs has been confirmed by 
Neuhauser et al (1999), with all the objects belonging to young clusters or star-forming 
regions. 

All these observations show that there is no drop in L x /L bo i in field stars up to a 
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spectral type of M7, which is well below the limit m ~ 0.3 M Q where stellar interiors are 
predicted to become fully convective. However, recent observations of the 2MASS objects 
(Kirkpatrick et al 1999, Fig. 15a) seem to show a significant decrease of activity for L-type 
objects, which confirms the decline of activity suggested previously for very late M-spectral 
types ( >M8)(Hawley et al 1996, Tinney & Reid 1998). In any event, X-ray emission does 
not disappear in fully convective stars and even at least young SSOs can support magnetic 
activity. A fossil field can survive only over a timescale ~ R 2 /i] ~ a few years for fully 
convective stars (77 = rj t is the turbulent magnetic diffusivity and R the stellar radius) so 
that a dynamo process is necessary to generate the magnetic field. The data suggest that 
(i) the dynamo-generation believed to be at work in the Sun does not apply for very-low 
mass objects, (ii) dynamo generation in fully convective stars is as efficient as in stars 
with a radiative core and (iii) whatever the dynamo mechanism is, it is bounded by the 
saturation condition at least for late-type stars. 

It is generally admitted that magnetic activity in the Sun results from the generation 
of a large-scale toroidal field by the action of differential rotation on a poloidal field at the 
interface between the convective envelope and the radiative core, where differential rotation 
is strongest, known as the tachocline (Spiegel & Weiss 1980, Spiegel & Zahn 1992). In 
this region a sufficiently strong magnetic field is stable against buoyancy. In this case the 
shear dominates the helicity, this is the so-called afl dynamo-generation, which predicts a 
correlation between activity and rotation, as observed in solar-type stars. 

As shown initially by Parker (1975), buoyancy prevents the magnetic field generated by 
the turbulent motions in a convective zone to become global in character. For this reason, 
the absence of a radiative core, i.e. of a region of weak buoyancy and strong differential 
rotation, precludes in principle the generation of a large-scale magnetic field. A dynamo 
generated by a turbulent velocity field that would generate chaotic magnetic fields in the 
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absence of rotation has been proposed as an alternative process for fully convective stars 
(Durney et al 1993). They found that a certain level of magnetic activity can be maintained 
without the generation of a large-scale field. The turbulent velocity field can generate 
a self-maintained small-scale magnetic field providing the magnetic Reynolds number 
Re = {vl/rf) 2 is large enough. The scale of this field is comparable to that of turbulence 
and in rough energy equipartition with the fluid motions i.e. ^{(Bp + Bfy) ~ \{pv1 onv ). 
Rotation is not essential in this case but it increases the generation rate of the field. 

Recently Kiiker and Chabrier (in preparation) explored an other possibility, namely the 
generation of a large-scale field by a pure a 2 -effect. In the a 2 -dynamo, helicity is generated 
by the action of the Coriolis force on the convective motions in a rotating, stratified fluid. 
The a 2 -effect strongly depends on the Rossby number, or the equivalent Coriolis number 
Q* = 2t conv Q = Ati/Ro. In low-mass objects, the convective turnover time is longer than 
the rotation period (see §2.1.3), so that VL* » 1. They find that the a 2 -dynamo is clearly 
supercritical. It generates a large-scale non-axisymmetric steady (co-rotating) field that 
is symmetric with respect to the equatorial plane. Equipartition of energy yields field 
strengths of several kiloGauss. The possible decrease of activity in very late-type stars, 
as observed for example with the small chromospheric activity of the M9.5 BD candidate 
BRI 0021-0214 in spite of its fast rotation (v sini =40 km s" 1 ) (Tinney et al 1998), is 
a more delicate issue that requires the inclusion of dissipative processes. Indeed, these 
calculations show that a 2 -dynamo can efficiently generate a large-scale magnetic field in the 
interior of fully convective stars. Since conductivity decreases in the outermost layers of the 
star, however, there no current would be created by the field in these regions, and thus no 
dissipative process and no activity (see e.g. Meyer & Meyer- Hofmeister 1999). 

The observed continuous transition in rotation and activity at the fully convective 
boundary suggests in fact that the a 2 -dynamo is already at work in the convection zones 
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of the more massive stars. An interesting possibility for verifying the present theory would 
be Doppler imaging of fast-rotating LMS and BDs. Turbulent dynamo is likely to yield a 
spatially uniform chromospheric activity whereas the large-scale a 2 process suggested by 
Kiiker & Chabrier generates asymmetry. Moreover, non-axisymmetric fields can propagate 
in longitudinal directions without any cyclic variation of the total field energy whereas 
dynamo waves generated by aVt processes propagate only along the lines of constant 
rotation rate (Parker 1955). Therefore we do not expect cycles for uniformly rotating (fully 
convective) stars. 

3. MECHANICAL AND THERMAL PROPERTIES 

3.1. Mechanical Properties 

Figure |3| portrays the mass-radius behaviour of LMS and isolated SSOs from the 
Sun to Jupiter for t = 6 x 10 7 (dash-dot line) and 5xl0 9 yr (solid line) for Z = Z & and 
Z = 10~ 2 x Z & (dash-line). The time required to reach the ZAMS, arbitrarily defined as 
the time when L nuc = 95% L tot , for solar metallicity LMS is given in Table III 

The general m-R behaviour reflects the physical properties characteristic of the 
interior of these objects, as inferred from Figure [I]. For m>O.2M , ip > 1 for all ages 
so that the internal pressure is dominated by the classical perfect gas ion+electron 
contribution (P = pkT / finiH, neglecting the molecular rotational/vibrational excited 
level contributions) plus the correcting Debye contribution arising from ion and electron 
interactions (Pdh — p 3 / 2 /T 1 / 2 ) and R oc m in first order from hydrostatic equilibrium. 
When the density becomes large enough during pre-MS contraction so that ip < I, which 
occurs at t > 50 Myr for m = 0.15 M and at t > 10 Myr for the HBMM m = 0.075 M , 
the EOS starts to be dominated by the contribution of the degenerate electron gas 
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(P oc p 5 / 3 ). This yields a minimum in the m-R relationship R mm ~ 0.08i? Q for m ~ 
0.06-0.07 M & . Full degeneracy (ip ~ 0) would yield the well-known zero-temperature 
relationship R oc m -1 / 3 , as in white dwarf interiors, but partial degeneracy and the 
non-negligible contribution arising from the (classical) ionic Coulomb pressure (which 
implies R oc m +1 / 3 ) combine to yield a smoother relation R = Rqui" 1 ' 8 at t=5 Gyr, 
where Rq ~ 0.06 R@ for 0.01 M Q < m< 0.07 M & , i.e. an almost constant radius. For the 
age and metallicity of the solar system this radius is of the order of the Jupiter radius 
Rj ~ 0.10 x R Q . At 5 Gyr, the radius reaches a maximum R ~ 0.11 R® for m~4x Mj 
(where Mj = 9.5 x 10 -4 M is the Jupiter mass) (Zapolsky & Salpeter 1969, Hubbard 1994, 
Saumon et al. 1996). Below this limit, degeneracy saturates (see Figure [I]) and the classical 
ionic pressure contribution becomes important enough so that we recover a nearly classical 
behaviour. Figure |^ also displays the astrophysically determined radii of the two eclipsing 
binary systems YY-Gem (Leung & Schneider 1978) and CM-Dra (Metcalfe et al 1996), 
and of the white dwarf companion GD448 (Maxted et al 1998). The bump on the 60 Myr 
isochrone near m ~ O.O1M results from initial deuterium burning (see §3.2). 

The m-R relationship is essentially determined by the EOS and is weakly sensitive 
to the outer boundary condition and the structure of the atmosphere, since the latter 
represents at most (for 1 M ) a few percent of the total radius of LMSs and SSOs (Dorman 
et al 1989, Chabrier & Baraffe 1997). The effect of metallicity Z on the m— R relation 
remains modest. A decrease in metallicity yields a slight decrease in the radius at a 
given age or at the same stage of nuclear burning (i.e the same H content). Indeed lower 
metallicity yields a larger T e $ at a given mass (see §3.2). This in turns implies an increase 
in the nuclear energy production to reach thermal equilibrium, and thus a larger central 
temperature T c . Since R oc fim/T from hydrostatic equilibrium, where /i is the mean 
molecular weight, determined essentially by hydrogen and helium, a lower-metallicity star 
contracts more to reach thermal equilibrium. As shown by Beuermann et al (1998), the 
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variation of R with Z predicted by the models of Baraffe et al (1998) is in agreement with 
the radii deduced from observations. 

3.2. Thermal Properties 

Figure |] exhibits the mass-effective temperature relationships for representative ages 
and metallicities. The arrows indicate the onset of formation of H 2 near the photosphere 
(Auman 1969, Copeland et al 1970, Kroupa et al 1990). As discussed in §2.2.2, molecular 
recombination favors convective instability in the atmosphere. Convection yields a smaller 
T-gradient (V ~ V a ^) and thus a cooler structure in the deep atmosphere (Brett 1995, 
Allard & Hauschildt 1995, Chabrier & Baraffe 1997 Figure 5). Therefore, a model with 
atmospheric convection corresponds to a larger T e g- since the (P, T) interior-atmosphere 
boundary is fixed for a given mass (§2.2.2). Furthermore, the adiabatic gradient in the 
regions of H 2 recombination decreases to a minimum value V a d ~ 0.1, compared to 
V a d ~ 0.4 for an ideal monoatomic gas (Copeland et al 1970, Saumon et al 1995 Figure 
17). Therefore, even if the atmosphere were already convective without H2, molecular 
recombination yields a flatter inner temperature gradient in the atmosphere, and thus 
enhances the former effect, i.e. a larger T e g for a given mass. Formation of H 2 in the 
atmosphere occurs at higher T e g for decreasing metallicity, because of the denser (more 
transparent) atmosphere (see Figure f|). The relation between the central and effective 
temperatures for LMS and SSOs can be inferred from Figures [l] and |] and is commented 
on by Chabrier and Baraffe (1997 §4.2) and Baraffe et al. (1998 §2). As these authors 
demonstrate, a grey approximation or the use of a radiative T(r) relationship significantly 
overestimates the effective temperature (from ~ 50 to 300 K depending on the atmospheric 
treatment). 

A representative T c -p c diagram is shown in Figure |5] (see also Figure 8 of Burrows et 
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al 1997) and illustrates the different evolutionary paths for LMS and SSOs in this diagram. 
The bumps appearing on the isochrones between 10 6 and 10 s yr and log T c ~ 5.4-5.8 result 
from initial deuterium burning. For stars, T c and p c always increase with time until they 
reach the ZAMS. For BDs, T c first increases for ~ 10 7 -10 9 yr, for masses between ~ 0.01 
and ~ 0.07 M & , respectively. Then, when degeneracy becomes dominant, T c reaches a 
maximum and decreases. For objects with m < 5 x 10~ 3 M Q , T c always decreases for t > 1 
Myr. In the substellar regime there is no steady nuclear energy generation, by definition, 
and the evolution is governed by the change of internal energy J M ^§dm and the release 
of contraction work J M ■j^&dm. Even when degeneracy becomes important, SSOs keep 
contracting, though very slowly, within a Hubble time. 

4. EVOLUTION 

4.1. Evolutionary tracks 

Figure ||] exhibits T e g(t) obtained from consistent non-grey calculations for several 
masses, for [M/H]fl=0 (helium mass fraction Y = 0.275) and [M/H] = -2.0 (Y = 0.25), 
respectively. Initial deuterium burning (with an initial mass fraction [D] = 2 x 10~ 5 ) for 
masses m> 0.013 M proceeds quickly, at the early stages of evolution, and lasts about 
~ 10 6 -10 8 years. Objects below this limit are not hot enough to fuse deuterium in their core 
(see Figure §). For masses above 0.07 M for [M/H]=Q and 0.08 M Q for [M/H]=-2.0, the 
internal energy provided by nuclear burning quickly balances the gravitational contraction 
energy, and after a few Gyr the lowest-mass star reaches complete thermal equilibrium 
(L — J e nuc dm, where e nuc is the nuclear energy rate), for both metallicities. The lowest 
mass for which thermal equilibrium is reached defines the HBMM and the related 
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hydrogen-burning minimum luminosity HBML. Stars with m > O.4M develop a convective 
core near the ZAMS for a relatively short time, depending on the mass and metallicity, 
which results in the bumps for 0.6 and 1 M at respectively ~ 10 8 yr and ~ 3 x 10 7 yr 
(Chabrier and Baraffe 1997, §3.2). The dotted lines portray T e g(t) for objects with solar 
abundance when grain opacity is taken into account in the atmosphere (see §4.5.2). The Mie 
opacity due to the formation of refractory silicate grains produces a blanketing effect that 
lowers the effective temperature and luminosity at the edge of the main sequence, an effect 
first noticed by Lunine et al (1989). However as a whole, grain formation only moderately 
affects the evolution near and below the bottom of the main sequence and thus the HBMM 
and HBML. Models with grainless atmospheres yield m = 0.072 M , L = 5 x 1O _5 L and 
T e fj = 1700 K at the H-burning limit, whereas models with grain opacity give m = O.O7M , 
L = 4 x 10~ 5 L & and T e fj = 1600 K, for solar composition (Chabrier et al 2000). As shown in 
Figure |6], young and massive BDs can have the same effective temperature (or luminosity) 
as older very-low-mass stars, a possible source of contamination for the determination of 
the local stellar luminosity function at the bottom of the MS (see §5). 

Note the quick decrease of T e g (and L) with time for objects below the HBMM, L oc t a 
with a ~ —5/4 (Stevenson 1991, Burrows et al 1994, 1997), with a small dependence of 
a on the presence of grains. Slightly below ~ 0.072 M (resp. 0.083 M ) for [M/H]=0 
(resp. [M/H] < -1), nuclear ignition still takes place in the central part of the star, but 
cannot balance the ongoing gravitational contraction (see Figure |2|). Although these objects 
are sometimes called "transition objects" we prefer to consider them as massive BDs, 
because strictly speaking they will never reach thermal equilibrium. Indeed, just below the 
HBMM, the contributions from the nuclear energy source / e nVLC dm and the entropy source 
/ T^-dm are comparable and cooling proceeds at a much slower rate than mentioned above. 
Below about O.O7M (resp. 0.08 M ) for [M/H]=0 (resp. [M/H] < -0.5), the energetic 
contribution arising from hydrogen-burning, though still present for the most massive 
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objects, is orders of magnitude smaller than the internal heat content, which provides 
essentially all the energy of the star (e nuc << T|^|). Once on the ZAMS the radius for 
stars is essentially constant, whereas for BDs the contraction slows down when ip <0.1 
(Figures [l] and ||). 

The effects of metallicity on the atmosphere structure (Brett 1995, Allard and 
Hauschildt 1995) and on the evolution (Saumon et al 1994, Chabrier and Baraffe 1997) 
have been discussed extensively in the previously-cited references and can be apprehended 
with simple arguments: the lower the metallicity Z, the lower the mean opacity R and 
the more transparent the atmosphere so that the same optical depth lies at deeper levels 
and thus higher pressure (4p = |). Therefore for a given mass (log g) the (T,P) interior 
profile matches, for a given optical depth r, an atmosphere profile with larger T e g (dash-line 
on Figure |]), and thus higher luminosity L since the radius depends only weakly on the 
metallicity. The consequence is a larger HBMM for lower Z since a larger L requires more 
efficient nuclear burning to reach thermal equilibrium, and thus a larger mass. 

A standard way to display evolutionary properties is a theoretical Hertszprung-Russell 
diagram (HRD). Such a HRD for LMS and SSOs is shown in Figure [F]for several masses and 
isochrones from 1 Myr to 5 Gyr, whereas Figure [8] shows evolutionary tracks in a log g-T e s 
diagram (see also Burrows et al 1997). These figures allow the determination of the mass 
and age of an object from the gravity and effective temperature inferred from its spectrum. 



4.2. Mass-magnitude 

One of the ultimate goals of stellar theory is an accurate determination of the mass of 
an object for a given magnitude and/or color. Figure |^ (see also Baraffe et al 1998) shows 
the comparison between theory and observationally-determined masses in the -fT-band. The 
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solid line corresponds to a 5xl0 9 yr isochrone, for which the lowest-mass stars have settled 
on the MS, for a solar metal-abundance ([M/H] = 0), whereas the dotted line corresponds 
to a 10 8 yr isochrone for the same metallicity and the dashed line corresponds to a 10 10 yr 
isochrone for [M/H] = —0.5, representative of the thick-disk population. A striking feature 
is the weak metallicity-dependence in the K-band, compared to the strong dependence in 
the V-band (see Figure 3 of Baraffe et al 1998). As discussed in this paper, this stems from 
two different effects. On one hand, the increasing opacity in the optical, dominated by TiO 
and VO lines, and the decreasing H 2 opacity in the K-band with increasing metallicity shift 
the peak of the flux toward larger wavelengths. Thus, for fixed T cS the V-flux decreases 
and the K-flux increases with increasing [M/H]. On the other hand, for a given mass, the 
total flux (and T cS ) decreases with increasing metallicity, as mentioned in §4.1. These two 
effects add up in the V-band and yield an important variation of the V-flux with metallicity. 
In the K-band, they cancel and yield similar fluxes for a given mass below ~ 0.4 M Q 
(T e ff ~ 3500-4300 K, depending on [M/H]) at different metallicities. These arguments 
remain valid as long as H 2 CIA does not significantly depress the K-band, as it does for 
very metal-depleted objects. 



4.3. Mass-Spectral type 

The knowledge of spectral type (Sp) is extremely useful for analyzing objects with 
unknown distance or with colors altered by reddening, as in young clusters, for example. 
The determination of a mass-Sp relationship provides a powerful complement to the 
mass-magnitude relationship for assigning a mass to an object. Based on spectroscopic 
observations of low-mass nearby composite systems, Kirkpatrick & McCarthy (1994) have 
determined an empirical m-Sp relationship for M-dwarfs, restricted to M2-M6 spectral 
types. A theoretical relation has been derived by Baraffe and Chabrier (1996) from M0 to 
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M10, emphasizing the non-linear behavior of this relationship near the bottom of the MS, 
for >M6. Figure [K] shows the m-Sp relationship for M-, K- and G-dwarfs based on the 
Baraffe et al (1998) models. The spectral type for K- and G-dwarfs is derived from the 
model synthetic color (J — K) through the empirical Sp(I — K) relation recently derived by 
Beuermann et al (1998). As stressed by Baraffe and Chabrier (1996), the m-Sp relationship 
depends crucially on age and metallicity. At ~ 1 Gyr, objects with Sp later than M10 are 
below the HBMM and can be considered as bona-fide BDs. This limit decreases to ~ M7 
at 100 Myr and ~ M6 at 10 Myr. A difficulty arises when dealing with very young objects 
(t < 10 Myr) because of gravity effects. No reliable Sp-color relationship is presently 
available for these objects for which spectroscopic and photometric properties are found to 
be intermediate between giants and dwarfs (Luhman 1999). The present m-Sp relationships 
should be extended in the future to dwarfs cooler than M9.5-M10, the so-called L-dwarfs 
and methane-dwarfs. Note however that SSOs evolve at different rates through a series of 
spectral types as they cool, so we cannot associate a given spectral type with a specific 
mass for these objects. 



4.4. Irradiated planets 

The mass-radius relation and evolutionary sequences described above reflect the 
relations for isolated objects. After a rapid initial accretion phase and subsequent 
hydrodynamical collapse, planets orbiting stars will evolve differently. As shown by 
Hubbard (1977), illumination from a parent star will yield thermal expansion of the less 
massive (low gravity) objects, toward an asymptotic temperature T eq set only by the 
thermalized photons from the parent star and toward a (larger) asymptotic equilibrium 
radius R p , with R p = 2a(L e9 /L*) 1 / 2 [l/(l - A)} 1/2 and T eq = (1 - A) 1/4 (-R* /2a) 1/2 T ef fir , 
where a is the orbital distance, A the Bond albedo, L* = AnaRlTef^ the parent star 
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luminosity, effective temperature and radius and L eq = AiraRpT^ denotes the equilibrium 
luminosity (see Saumon et al 1994). As noted in Guillot et al (1996) and Guillot (1999), 
EGPs in close orbits are heated substantially by their parent star and their atmosphere 
cannot cool substantially. They develop an inner radiative region, as a result of this stellar 
heating and contract at almost constant T e Q, at a much smaller rate than if they were not 
heated by the star. Planets whose luminosity is larger than the absorbed stellar flux of the 
parent star just evolve along a fully convective Hayashi track for ~ 10 6 yr before reaching 
the afore-mentioned equilibrium state. 

The flux from an irradiated planet includes two contributions: the intrinsic thermal 
emission and the reflected starlight (the Albedo contribution) (see e.g. Seager & Sasselov 
1998, Marley et al 1999): 



k = (^) 2 ^ + (j) pm^fi^fK (5) 

a 4 da 

where d is the distance of the system to Earth and P(4>) is the dependence of the reflected 
light upon the phase angle between the star, the planet and Earth (P = 1 if the light 



reflected by the planet is redistributed uniformly over An steradians). Figure |TT| illustrates 
the flux from a young EGP (assuming the irradiation is isotropic) with T e g- = 1000 K 
orbiting a Sun-like star at 15 pc for various orbital distances. 



4.5. Color-magnitude diagrams 

Thanks to the recent progress both on the observational side, with the development of 
several ground-based and space-based observational surveys of unprecedented sensitivity in 
various optical and infrared passbands, and on the theoretical side, with the derivation of 
synthetic spectra and consistent evolutionary calculations based on non-grey atmosphere 
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models, it is now possible to confront theory with observation directly in the observational 
planes, i.e. in color-color and color-magnitude diagrams (CMD). This avoids dubious 
color-T fj or color- Mboi transformations and allows more accurate determinations of the 
intrinsic properties (m, L, T e ff, R) of an object from its observed magnitude and/or color. In 
this section we examine the behaviour of LMS and SSOs in various CMDs characteristic of 
different populations in terms of age and metallicity. These diagrams capture the essence of 
the observational signatures of the very mechanical and thermal properties of these objects, 
described in the previous sections. 

4-5.1. Pre-main sequence and young clusters 

Numerous surveys devoted to the detection of SSOs have been conducted in young 
clusters with ages spanning from ~ 1-10 Myr to > 10 8 yr for the Pleiades or the Hyades 
(see Martin 1999 and Basri this volume for reviews). Observations of young clusters present 
two important advantages, namely (i) all objects in the cluster are likely to be coeval 
within a reasonable range, except possibly in star forming regions, where the spread in 
ages for cluster members can be comparable with the age of the cluster, (ii) young objects 
are brighter for a given mass (see Figure |6|), which makes the detection of very-low mass 
objects easier. Conversely, they present four major difficulties, (i) extinction caused by the 
surrounding dust modifies both the intrinsic magnitude and the colors of the object, (ii) 
accurate proper motion measurements are necessary to assess whether the object belongs 
to the cluster, (iii) gravity affects both the spectrum and the evolution, (iv) the evolution 
and spectrum of very young objects (t< 1 Myr) may still be affected by the presence of 
an accretion disk or circumstellar material residual from the protostellar stage. Young 
multiple systems remove some of these difficulties and provide excellent tests for models 
at young ages. Like for example the quadruple system GG TAU (White et al 1999), with 
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components covering the whole mass-range of LMS and BDs from 1 M Q to ~ 0.02 M Q . 
Models based on non-grey atmospheres (Baraffe et al 1998) are the only ones consistent 
with these observations (White et al 1999, Luhman 1999). The comparison of ~ 10 6 yr 
isochrones for SSOs with observations in a near-IR CMD is shown in Zapatero Osorio et al 
(1999) for the young cluster a-Orionis (their figure 1). 



4-5.2. Disk field stars 



Figures [12] and [13] display the CMDs for LMS and SSOs in the optical and the infrared 
for various ages and metallicities (see also Burrows et al 1997). Various sets of models 
correspond to calculations where (i) grain formation is included in the atmosphere EOS but 
not in the transfer equations, mimicking a rapid settling below the photosphere (COND 
models), (ii) grain formation is included both in the EOS and the opacity (DUSTY models) 
(see Allard et al 1999), (iii) grainless models. Models (i) and (ii) represent extreme cases 
that bracket the complex grain processes at play in these atmospheres (see §2.2.1). The 
optical sequence shows a monotonic magnitude-color behaviour with 3 changes of slope 
around My ~ 10, 13 and 19 for Z = Z Q , respectively, which correspond to 0.5 M , 
T eff =3600 K; 0.2 M , T eff =3300 K and 0.08 M Q , T cS = 2330 K on the 5 Gyr isochrone. 
They reflect, respectively, the onset of convection in the atmosphere, degeneracy in the core 
and grain absorption near the photosphere, as described in the previous sections. The latter 
feature, which yields a steep increase of My at (V — I) ~ 5 does not appear in dust-free 
models (cf. Baraffe et al 1998), where TiO is not depleted by grain formation in the 
atmosphere and thus absorbs the flux in the optical. Interestingly enough, this steepening is 
observed in the sample collected recently by Dahn et al (1999) which includes some DENIS 
and 2MASS L-dwarf parallaxes. This effect is also observed in the (R — /) CMD of the 
Pleiades (Bouvier et al 1998), in agreement with the theoretical predictions (Chabrier et al 
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2000). This clearly illustrates the effect of grain formation in the atmosphere which affects 
the spectra of LMS and SSOs below T e g- < 2400 K. In the near-infrared color (J-K) (Figure 



13), grain condensation yields very red colors for the DUSTY models, whereas COND and 
grainless models loop back to short wavelengths (blueward) below m ~ 0.1LM©, Mr- ~ 10 
for [M/H] = -2, and m ~ 0.072M o (resp. ~ 0.04M Q ) at 5 Gyr (resp. 100 Myr) for solar 
metallicity, i.e. M# ~ 12. As described in §2.2, this stems from (i) the onset of H 2 CIA 
for low-metallicity objects and (ii) the formation of CH4 at the expense of CO when the 
peak of the Planck function moves in the wavelength range characteristic of the absorption 
bands of this species, for solar-like abundances. 

In terms of colors there is a competing effect between grain and molecular opacity 
sources for objects at the bottom and just below the MS. The backwarming effect resulting 
from the large grain opacity destroys molecules such as e.g H 2 0, one of the main sources 
of absorption in the near-IR (cf. §2.2.1), and yields the severe reddening near the bottom 
of the MS (T eff <2300 - 2400 K), as illustrated by the DENIS (Tinney et al 1998) and 
2MASS (Kirkpatrick et al 1999a) objects and by GD165B (Kirkpatrick et al 1999b). As the 
temperature decreases, grains condense and settle at deeper layers below the photosphere 
while methane absorption in the IR increases (see §2.2) so that the peak of the flux will 
move back to shorter wavelengths and the color sequence will go from the DUSTY one to 
the COND one. The exact temperature at which this occurs is still uncertain because it 
involves the complex grain thermochemistry and dynamics outlined in §2.2, but it should 
lie between T cS ~ 1000 K, for G1229B, and ~ 1700 K, corresponding to J-K ~ 2 for the 
presently reddest observed L-dwarf with determined parallax LHS102B (Goldman et al 
1999, Basri et al 2000). 

In spite of the afore-mentioned strong absorption in the IR, BDs around 1500 K radiate 
nearly 90% (99% with dust) of their energy at wavelengths longward of 1 /im and infrared 
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colors are still preferred to optical colors (at least for solar metal abundance), with J, Z, H 
as the favored bands for detection, and Mm ~M' l ~ 10 — 11, ~Mj ~ 11 — 12 at the 
H-burning limit, at 1 Gyr (see Chabrier et al 2000). 

4-5.3. Halo stars. Globular clusters 

Observations of LMS belonging to the stellar halo population of our Galaxy (also 
called "spheroid" to differentiate it from the p(r) oc 1/r 2 dark halo) or to globular clusters 
down to the bottom of the MS have been rendered possible with the HST optical (WFPC2) 
and IR (NICMOS) cameras and with parallax surveys at very faint magnitudes (Leggett 
1992, Monet et al 1992, Dahn et al 1995). This provides stringent constraints for our 
understanding of LMS structure and evolution for metal-depleted objects. The reliability of 
the LMS theory outlined in §2 has been assessed by the successful confrontation to various 
observed sequences of globular clusters ranging from [M/H] ~ -2.0 to -1.00 (Baraffe et al 
1997). Figure 8 of Baraffe et al (1997) portrays several tracks for different metallicities, 
corresponding to the MS of various globular clusters. The tracks are superposed with the 
subdwarf sequence of Monet et al (1992), which is identified as stellar halo objects from 
their kinematic properties. As shown in this figure, the pronounced variations of the slope 
in the CMD around ~ 0.5 M & and ~ 0.2 M , which stem from the very physical properties 
of the stars, namely the onset of molecular absorption in the atmosphere and of degeneracy 
in the core (see §2), are well reproduced by the theory at the correct magnitudes and colors. 
The predicted blue loop in IR colors caused by the ongoing CIA absorption of H 2 for the 

2 Globular clusters with an observed [Fe/H] must be compared with theoretical models 
with the corresponding metallicity [M/H] = [Fe/H] + [O/Fe], in order to take into account 
the enrichment of a-elements (see Baraffe et al 1997) 
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lowest-mass (coolest) objects (see §2.2.1 and Figure [X3| for [M/H] = —2), has also been 
confirmed by observations with the NICMOS camera (Pulone et al 1998). 

4.6. The lithium and deuterium tests 

One ironclad certification of a BD is a demonstration that hydrogen fusion has not 
occurred in its core. As Figure |2| illustrates, lithium burning through the Li 7 (p, a)He 4 
reaction occurs at a lower temperature than is required for hydrogen fusion. The timescale 
for the destruction of lithium in the lowest-mass stars is about 100 Myr. Moreover the 
evolutionary timescale in LMS /BDs is many orders of magnitude larger than the convective 
timescale (see e.g. Bildsten et al 1997), so core abundances and atmospheric abundances 
can be assumed identical. Therefore, the retention of lithium in a fully mixed object older 
than 10 8 yr signifies the lack of hydrogen burning. This provides the basis of the so-called 
"lithium-test" first proposed by Rebolo et al (1992) to identify bona-fide BDs. Evolutionary 
models based on non-grey atmospheres and the screening factors described in §2.1.2 yield a 
lithium-burning minimum mass ran ~ 0.06 M & (Chabrier et al 1996). Note from Figure |2] 
the strong age dependence of the lithium-test : young stars at t < 10 8 yr (depending on the 
mass) will exhibit lithium (which in passing precludes the application of the lithium-test for 
the identification of SSOs in star forming regions), whereas massive BDs within the mass 
range [0.06-0.07 M & ] older than ~ 10 8 yr will have burned lithium. 

By observing a young cluster older than ~ 10 s yr, one can look for the boundary- 
luminosity below which lithium has not yet been depleted. Objects fainter than this limit 
will definitely be in the substellar domain. Conversely, the determination of the lithium 
depletion edge yields the age of the cluster, as was first proposed by Basri et al (1996) (see 
also e.g. Stauffer et al 1998). Indeed, the lithium test provides a "nuclear" age that may 
be even more powerful a diagnostic than the conventional nuclear age from the upper main 
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sequence turnoff, since the evolution of the present fully convective low-mass objects does 
not depend on ill-constrained parameters such as mixing length or overshooting. As noted 
by Stauffer et al (1998) the age-scale for open clusters based on the (more reliable) lithium 
depletion boundary has important implications for stellar evolution: the ages of several 
clusters (Pleiades, a-Per, IC 2391) are all consistent with a small but non-zero amount of 
overshooting to be included in the evolutionary models at the turn-off mass in order to 
yield similar ages. See Basri (this issue) for a more detailed discussion of the lithium-test. 

The deuterium-test can be used in a similar manner, extending the lithium test to 
smaller masses and younger ages, typically t<10 7 yr (Bejar et al 1999). The D-burning 
minimum mass is predicted as mc ~ 0.013M Q (see Figure 0) and has been proposed 
(Shu et al 1987) as playing a key role in the formation of isolated star-like objects, in 
contrast to objects formed in a protoplanetary disk. Deuterium-depletion is illustrated in 
Figure [14]. This figure displays the evolution of SSOs above mp in the I-band. The left 
and right diagonal solid lines correspond to 50% D-depletion ([£)]/[£)] = 1/2) and 99% 
D-depletion, respectively ([D]q = 2 x 10 -5 is the initial mass fraction). The inset displays 
the corresponding curves in T eff . Spectroscopic signatures of deuterium include absorption 
lines of deuterated water HDO between 1.2 and 2.1 jum (Toth 1997). The observation of 
deuterated methane CH 3 D can not be used for the deuterium test clS clll cl ee-indicator: at 
the temperature methane forms ( < 1800 K), SSOs above the D-burning minimum mass are 
old enough for all deuterium to have been burned (see Figure . 



4.7. Low- mass objects in cataclysmic variables 

The study of cataclysmic variables (CVs) is closely related to the theory of LMS 
and SSOs regarding their inner structure and the magnetic field generation. CVs are 
semi-detached binaries with a white dwarf (WD) primary and a low-mass stellar or 
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sub-stellar companion (m < 1 M©) that transfers mass to the WD through Roche-lobe 
overflow (see e.g. King 1988 for a review). In the current standard model, mass transfer 
is driven by angular momentum loss caused by magnetic wind braking for predominantly 
radiative stars and by gravitational wave emission for fully convective objects. The Roche 
lobe filling secondary's mean density p determines almost entirely the orbital period P, with 
Ph — k/p 1 / 2 (where k ~ 8.85 is a weak function of the mass ratio, Ph is the orbital period 
in hr, and p = (mz/ R\) / {M & / R Q ) is the mean density of the secondary in solar units) (see 
King 1988). For objects obeying classical mechanics, i.e. in the stellar regime, m oc R (see 
§3.1 and Figure |3|) so that the mean density increases with decreasing mass (p oc 1/m 2 ) and 
the orbital period decreases along evolution as the secondary loses mass Ph oc 1/p 1 ^ 2 oc m 2 . 
The situation reverses near and below the HBMM, when electron degeneracy dominates. 
In that case, the radius increases very slightly with decreasing mass (Figure ||), so that 
the mean density is essentially proportional to the mass (p oc m) and the orbital period 

— 1/2 

increases along evolution Ph oc~ m 2 . Then the secular evolution of the orbital period 
reverses when the donor becomes a BD. The analysis of CV orbital evolution thus provides 
important constraints on LMS and BD internal structure. 

The major puzzle of CVs is the observed distribution of orbital periods, namely the 
dearth of systems in the 2-3 h period range and the minimum period P m i n at 80 min (see 
King 1988). The most popular explanation for the period gap is the disrupted magnetic 
braking scenario (Rappaport et al 1983, Spruit & Ritter 1983), and most evolutionary 
models including this process do reproduce the observed period distribution. However, the 
recent progress realized in the field of LMS and SSOs now allows the confrontation of the 
theoretical predictions with the observed atmospheric properties (colors, spectral types) of 
the secondaries (Beuermann et al 1998, Clemens et al 1998, Kolb & Baraffe 1999a). The 
afore-mentioned standard period-gap model may be in conflict with the observed spectral 
type of some CV secondaries (Beuermann et al 1998, Kolb & Baraffe 1999b). Although 
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alternative explanations exist for the period gap, none of them has been proved as successful 
as the disrupted magnetic braking scenario. The most recent alternative suggestion, based 
on a characteristic feature of the mass-radius relationship of LMS (Clemens et al 1998), 
has been shown to fail to reproduce the observed period distribution around the period 
gap (Kolb et al 1998). On the other hand, as discussed in §2.3, no change in the level of 
activity is observed in isolated M-dwarfs along the transition between partially and fully 
convective structures at m ~ O.35M and Sp ~ M2-M4, which indicates that a magnetic 
field is still generated in the mass-range corresponding to the CV-gap. In this case, the 
abrupt decrease of angular momentum losses by magnetic braking at the upper edge of 
the period gap could result from a rearrangement of the magnetic field when stars become 
fully convective, without necessarily implying a sudden decline in the magnetic activity, as 
suggested by Taam & Spruit (1989, see also Spruit 1994). 

Finally, a ~ 10% discrepancy still remains between the theoretical and the observed 
minimum period P m i n =80 min, even when including the most recent improvements in 
stellar physics (see Kolb & Baraffe 1999a for details). Residual shortcomings in the theory, 
either in the EOS or in the atmosphere, cannot be ruled out as the cause of this discrepancy. 
Alternatively, an additional driving mechanism to gravitational radiation in fully convective 
objects can reconcile predicted and observed P min (Kolb & Baraffe 1999a). Because 
magnetic activity is still observed in fully convective late spectral type M-dwarfs, magnetic 
braking could operate in CV secondaries even down to P m i n , but with a weaker efficiency 
than above the period gap. These open questions certainly require a better understanding 
of magnetic field generation and dissipation in LMS and BDs, a point already stressed in 
§2.3. 



5. GALACTIC IMPLICATIONS 
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5.1. Stellar Luminosity Function and Mass Function 

It is now well established that visible stars are not numerous enough to account 
for the dynamics of our Galaxy, the so-called galactic dark matter problem. A precise 
determination of their density requires the correct knowledge of the luminosity function 
(LF) down to the H-burning limit, and a correct transformation into a mass-function (MF). 
The latter issue has been improved significantly with models that describe more accurately 
the color-color and color-magnitude diagrams of LMS and BDs for various metallicities, 
and most importantly that provide mass- magnitude relationships in good agreement with 
the observations (see §4). The former issue, however, is not completely settled at present 
and significant differences still exist between various determinations of the nearby LF and 
between the nearby LF and the photometric LF determined from the ground and with the 
HST, as shown later in this section. 

The main problem with the nearby sample is the limited number of stars at faint 
magnitudes. Although the sample is complete to 20 pc for stars with My < 9.5, it is severely 
incomplete beyond 5 pc for My > 12 (Henry et al 1997). Kroupa (unpublished) derived a 
nearby LF $ nea r by combining Hipparcos parallax data, which is essentialy complete for 
My < 12 at r=10 pc, and the sample of nearby stars (Dahn et al 1986) with ground-based 
parallaxes for V > 12 to a completeness distance r=5.2 pc. The nearby LF determined by 
Reid and Gizis (1997) is based on a volume sample within 8 pc. Most of the stars in this 
survey have parallaxes. For all the late K and M dwarfs, however, trigonometric parallaxes 
are not available and these authors use a spectroscopic TiO-index vs My l ° relation to 
estimate the distance (Reid et al 1995). This sample was revised recently with Hipparcos 
measurements and new binary detections in the solar neighborhood (Delfosse et al 1999a) 
and leads to a revised northern 8-pc catalogue and nearby LF (Reid et al 1999). These 
authors argue that their sample should be essentially complete for My < 14. However, the 
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analysis of completeness limits by Henry et al (1997, their figure 1) shows that the known 
stellar census becomes substantially incomplete for distances larger than 5 pc. About 
35% of the systems in the Reicl et al sample are multiple and ~ 45% of all stars have a 
companion in binary or multiple systems. 

The photometric LFs, & p hot, based on large observed volumes via deep pencil-beam 
surveys, avoid the limitation due to small statistics but introduce other problems such as 
Malmquist bias and unresolved binaries (see e.g. Kroupa et al 1993, Kroupa 1995), yielding 
only the determination of the stellar system LF. For the HST LF (Gould et al 1998), with 
/ < 24, the Malmquist bias is negligible because all stars down to ~ 0.1 M & are seen through 
to the edge of the thick disk. A major caveat of photometric LFs, however, is that the 
determination of the distance relies on a photometric determination from a color-magnitude 
diagram. In principle this requires the determination of the metallicity of the stars, since 



colors depend on metallicity (see Figure [12]). The illustration of the disagreement between 
the three afore-mentioned LFs, $5 pc , $8pc and &hst, at faint magnitude is apparent in 
Figure pl| The disagreement between the two nearby LFs, in particular, presently has 
no robust explanation. It might come from incompleteness of the 8 pc sample at faint 
magnitudes. Note also that in the spectroscopic relation used to estimate the distance, 
My can be uncertain by ~ 1 mag (see Figure 3 of Reid et al 1995), certainly a source 
of Malmquist bias, even though the number of stars without parallax in this sample is 
small (~ 10%) and a Malmquist bias on this part should not drastically affect the results. 
On the other hand, as noted in §4.1, the end of the 5-pc sample might be contaminated 
by a statistically significant number of young BDs, or by stars with slightly under-solar 
abundances. 



Figure 16 compares the MFs obtained from the two afore- mentioned nearby LFs, using 



the mass-magnitude relations described in §4.2. In practice, the observed sample constitutes 
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a mixture of ages and metallicities. However the Hipparcos CMD indicates that ~ 90% of 
disk stars have —0.3 < [M/H] < 0.1 (Reid 1999) so the spread of metallicity is unlikely to 
significantly affect the derivation of the MF through the mass-magnitude relation. Moreover 
the choice of m-M# minimizes metallicity effects (see Figure |9| and figure 3 of Baraffe et al 
1998). The MS lifetimes for stars with m< 1 M & are longer than a Hubble time and age 
variations in the mass-magnitude relation will affect only objects <O.15M younger than 
~ 0.5 Gyr (see Figure ^). 

Superposed is a power-law MF normalized at 0.8 M Q on the Hipparcos sample 
dn/dm = 0.031(m/0.8)~ a = 0.024 m~ a M~ 1 pc~ 3 with a = 1.2. As shown in Figure [IB], this 
provides a very reasonable description of the low-mass part of the MF, within a variation of 
~ 10% for a, thus confirming the previous analysis of Kroupa et al (1993), except for the 
very last bin obtained from Kroupa's LF, which predits about twice as many stars at 0.1 
M & (out of the figure). Similar results are obtained from the bolometric and Mk- LFs for 
the 8-pc sample. Integration of this MF yields very-low-mass stars (m < 0.8 M©) number- 
and mass- densities uvlms — 7.0 x 10~ 2 pc~ 3 and Pvlms ~ 2.0 x 10~ 2 M pc~ 3 , respectively. 
Adding up the contribution from more massive stars, p>o.s ~ 1.4 x 10 -2 M Q pc -3 (Miller & 
Scalo 1978) and stellar remnants, Pwd+ns ~ 3.0 x 10~ 3 M pc~ 3 yields a stellar density 
p± m 3.4 x 10 -2 M Q pc~ 3 in the Galactic disk, i.e. a surface density w 22 ± 2 M pc^ 2 
(assuming a scale height h = 320 pc). A Salpeter MF, dn/dm oc m -2 ' 35 , all the way down 
to the bottom of the MS would overestimate the density at 0.1 M Q by more than a factor 
of 10. 

For the Galactic spheroid the question is more settled. The LFs determined from 
nearby surveys (Dahn et al 1995) and from the HST (Gould et al 1998) are comparable 
within about 2 sigmas, and yield a MF with «<1 (Graff & Freese 1996, Chabrier & 
Mera 1997, Gould et al 1998) and a stellar number- and mass-density n± < 10 -3 pc~ 3 and 
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p± < 4.0 x M Q pc~ 3 , i.e. an optical depth r ~ 10~ 9 , about 1% of the value measured 
toward the LMC. At present, the dark halo MF is completely unknown but both the Hubble 
Deep Field observations and the narrow-range of the observed time-distribution of the 
microlensing events toward the LMC strongly suggest an IMF different from a Salpeter IMF 
below ~ IMq (Chabrier 1999). Assuming a homogeneous distribution, both constraints 
yield a mass-density for LMS in the dark halo about 2 orders of magnitude smaller than 
the afore-mentioned spheroid one (Chabrier & Mera 1997). This means that essentially no 
stars formed below about a solar mass in the dark halo. 

5.2. Brown dwarf mass function 

Even though the possibility that BDs could make up for the Galactic missing mass is 
now clearly excluded, a proper census of the number of BDs has significant implications 
for our understanding of how stars and planets form. The determination of the BD MF is 
a complicated task. By definition, BDs never reach thermal equilibrium and most of the 
BDs formed at the early stages of the Galaxy will have dimmed to very low-luminosities 
(L oc l/t). Thus observations will be biased toward young and massive BDs. This 
age-undetermination is circumvented when studying the BD MF in clusters because objects 
in this case are likely to be coeval. The Pleiades cluster has been extensively surveyed and 
several BDs have been identified down to ~ 0.04 M Q (Martin et al 1998, Bouvier et al 1998, 
Hambly et al 1999). A single power-law function from ~ 0.4 to 0.04 M seems to adequately 
reproduce the observations (not corrected for binaries) with some remaining uncertainties 
in the exponent: a ~ 0.6 — 1.0. However, stellar objects in very young clusters ( < 10 6 yr) 
might still be accreting whereas older clusters may have already experienced significant 
dynamical evolution, mass segregation in the core and/or evaporation in the outer regions 
(see e.g. Raboud & Mermilliod 1998). In this case the present-day mass function does not 
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reflect the initial mass function. For these reasons, it is probably premature to claim a 
robust determination of the MF in young clusters. 

The DENIS (Delfosse et al 1999b) and 2MASS (Kirkpatrick et al 1999b; Burgasser et 
al 1999) surveys, which have covered an area of several hundred square degrees and are 
complete to K — 13.5 and 14.5, respectively, have revealed about 20 field L-dwarfs and 4 
field methane-BDs. This yields an L-dwarf number- density til ~ 0.03 ± 0.01 sq deg -1 for 
K < 14.5 and a methane-BD number- density n C H 4 ~ 0.002 sq deg -1 for J < 16. Although 
these numbers correspond to small statistics and should be considered with caution, they 
provide the first observational constraints on the substellar MF in the Galactic disk and 
are consistent with a slowly rising BD MF with a ~ 1-2 (Reid et al 1999; Chabrier in 
preparation), assuming a constant formation rate. 

Independent, complementary information on the stellar and substellar MF comes 
from microlensing observations. Indeed, the time distribution of the events provides a 
(model-dependent) determination of the mass-distribution and thus of the minimum mass 
of the dark objects: dN ev /dt e = E x e(t e ) x dT/dt e oc P(m)/y / m, where E is the observed 
exposure, i.e. the number of star x years, e is the experimental efficiency, T is the event 
rate and P(m) is the mass probability distribution. The analysis of the published 40 
MACHO (Alcock et al 1997) + 9 OGLE (Udalsky et al 1994) events toward the bulge 
is consistent with a rising MF at the bottom of the MS whereas a decreasing MF below 
0.2 M Q seems to be excluded at a high (> 90%) confidence level (Han & Gould 1996, 
Mera et al 1998). Although the time distribution might be affected by various biases (e.g. 
blending) and robust conclusions must wait for larger statistics, the present results suggest 
that, in order to explain both star counts and the microlensing experiments, a substantial 
number of BDs must be present in the Galactic disk. Extrapolation of the stellar MF 
determined in §5.1 into the BD domain down to 0.01 M yields for the Galactic disk a BD 
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number-deiasitj comparable to the stellar one, ubd ~ 0.1 pc 3 ~ n+, and a mass-density 
Pbd ~ 3.0 x l(r 3 M pc- 3 , i.e. S BD « 2M pc~ 2 . 

For the spheroid, extrapolation of the previously-determined stellar MF yields 
fi BD < 10~ 4 pc~ 3 , pbd < lO~ 5 M pc~ 3 , less than 0.1% of the required dynamical density. 
For the dark halo the density is about 2 orders of magnitude smaller, as mentioned before. 

It is obviously premature to try to infer the mass distribution of exoplanets. This 
will first require a clear theoretical and observational distinction between planets and 
brown dwarfs. However, an interesting preliminary result comes from the observed mass 
distribution of the companions of G and K stars. As shown in figure 4 of Mayor et al 
(1998), there is a strong discontinuity in the mass distribution at m 2 /shu 5Mj, with a 
clear peak below this limit. This suggests that planet formation in a protoplanetary disk is 
a much more efficient mechanism than BD formation as star companions, at least around 
G and K stars (see e.g. Marcy & Butler 1998). It also suggests that the MF for BD stellar 
companions differs from the MF in the field. Ongoing observations around M-stars will tell 
us whether such a mass distribution still holds around low-mass stars. 



6. Conclusions 

This review has summarized the significant progress achieved within the past few 
years in the theory of cool and dense objects at the bottom of and beyond the main 
sequence: low-mass stars, brown dwarfs and gaseous planets. The successful confrontation 
of the theory with the numerous detections of low-mass stellar and sub-stellar objects 
allows a better understanding of their structural and thermal properties, and allows 
reliable predictions about their evolution. This in turn brings confidence in the predicted 
characteristic properties of these objects, a major issue in terms of search strategies 
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for future surveys. Important problems remain to be solved to improve the theory. A 
non-exhaustive list includes for example: (i) a better determination of the EOS in the 
pressure-ionization region, with the possibility of a first-order phase transition, (ii) the study 
of phase separation of elements in SSO interiors, (iii) an improved treatment of convection 
in optically-thin regions, (iv) a precise description of the dynamics of grain formation and 
sedimentation in SSO atmospheres, (v) the derivation of an accurate mass-T e g-age scale 
for SSOs and young objects, (vi) a correct understanding of magnetic-field generation and 
dissipation in active LMS and BDs. 

The increasing number of observed LMS and SSOs, together with the derivation of 
accurate models, eventually will allow a robust determination of the stellar and substellar 
mass functions, of the minimum mass for the formation of star-like objects, and thus of the 
exact density of these objects in the Galaxy. As we discussed in §5, present determinations 
in various Galactic regions point to a slowly rising MF near and below the H-burning 
limit, with a BD number- density comparable to the stellar one, and a MF truncated below 
~ 1 M© in the dark halo. A more precise determination must await confirmation from 
future observations. At last, the amazingly rapid pace of exoplanet discoveries should yield 
the determination of the planetary MF and maximum mass, and eventually the direct 
detection of such objects, a future formidable test for the theory. 
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Table 1: ZAMS age for solar metallicity LMS (after Baraffe et al (1998)) 
m/M Q 0.075 0.08 0.09 0.1 0.2 0.4 0.6 0.8 1 
tzAMS (Gyr) 3 1.7 1 0.7 0.33 0.20 0.12 0.06 0.04 
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Figure Legends 

Fig. 1. — Central temperature (in K), density (in g cm -3 ), and degeneracy parameter along 
the LMS-SSO mass range for objects with Z = Z Q at 5 Gyr (solid line) and 10 s yr (dashed 
line), and with a metallicity Z = 10~ 2 Z Q at 5 Gyr (dotted- line). 

Fig. 2. — Central temperature as a function of age for different masses. Th, Tu and Td 
indicate the hydrogen, lithium and deuterium burning temperatures, respectively. 

Fig. 3. — Mass-radius relationship for LMS and SSOs for 2 ages, t — 6 x 10 7 yr (dash-dot) 
and 5xl0 9 yr (solid) for Z = Z & , and t = 5 x 10 9 yr for Z = 1(T 2 x Z Q (dash). The 
HBMM is 0.075 M for Z = Z & and 0.083 M for Z = 10~ 2 x Z e . Also indicated are 
the observationally-determined radii of various objects (see text) and the position of Jupiter 
radius (J). The bump on the 6 x 10 7 yr isochrone illustrates the initial D-burning phase. 

Fig. 4. — Mass-T cfr relationship for LMS and SSOs. Same legend as in Figure 3. The arrows 
indicate the onset of H2 formation near the photosphere. Note the log scale in the inset. 

Fig. 5. — T c -p c relationship (in cgs) for LMS (solid lines) and SSOs (dashed lines) from 1 M 
to 0.001 M (masses in M indicated on the curves). Dotted lines represent 10 6 , 10 7 , 10 8 , 
10 9 and 5 10 9 yr isochrones from bottom to top. The bumps on the 10 6 -10 8 yr isochrones at 
log T c ~ 5.4 — 5.8 correspond to the initial deuterium burning phase. 

Fig. 6. — Effective temperature vs time (yr) for objects from 1 M to 10~ 3 M (masses 
indicated in M ). Solid lines: Z = Z Q , no dust opacity; dotted lines: Z = Z Q , dust opacity 
included, shown for 0.01, 0.04 and 0.07 M ; dashed line: Z = 10" 2 x Z Q (only for 0.3 M ). 
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Fig. 7. — Theoretical H-R diagram for various masses (labeled in M ). The weak dependence 
of radius on mass for SSOs yields the merging of the tracks for the lowest-mass objects. 
Dotted lines represent 10 6 , 10 7 , 10 s and 5 10 9 yr isochrones from right to left. The models 
are calculated with a mixing length I = 1.9 Hp. 

Fig. 8. — log g (in cgs) vs T eS (in K) for LMS (solid curves) and SSOs (dashed curves) from 
1 M G to 0.001 M Q (masses in M Q indicated on the curves). Dotted lines represent 10 6 , 10 7 , 
10 8 and 5 10 9 yr isochrones from bottom to top. 

Fig. 9. — Mass-M^ relationship. Observationally determined masses are indicated by filled 
triangles (Henry & McCarthy 1993) and circles (Delfosse et al 1999c, Forveille et al 1999). 
Solid line: [M/H] = 0, t = 5 x 10 9 yr; dotted line: [M/H] = 0, t = 10 8 yr; dash-line: 
[M/H] = -0.5, t = 10 10 yr. Note the log-scale in the inset. 



Fig. 10. — Mass-Spectral type relationships. Observations are from: Kirkpatrick & 
McCarthy (1994) for disk M-dwarfs (circles) and Clausen et al. (1999) for K- and G- dwarfs 
(squares). Note that the latter objects are likely to be young. The Sun is included at Sp = 
G2. The theoretical models are from Baraffe et al (1998). Solid line: [M/H] = 0, t = 5 x 10 9 
yr; dotted line: [M/H] = 0, t = 10 8 yr; dot-dash-line: [M/H] = 0, t = 3 x 10 7 yr; dash-line: 
[M/H] = -0.5, t = 5 x 10 9 yr. 
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Fig. 11. — Spectrum of a young EGP (T cff = 1000K, defined as the effective temperature 
of the non-irradiated object) irradiated by a G2V primary at 0.2,0.1, and 0.05 AU distance 
(bottom to top). The plot shows the incident spectrum (topmost curve), the spectrum 
emitted by the irradiated planet (magenta curves with more flux and shallower water and 
methane bands) and the non-irradiated spectrum of the EGP (bottom full black curve) 
assuming a distance of 15 pc to Earth. The dashed black line is the blackbody spectrum 
for T cff = 1000K. The EGP model includes dust opacities, (courtesy of F. Allard and P.H. 
Hauschildt) 

Fig. 12. — My vs (V-I) diagram for different ages and metallicites: DUSTY models (see 
text) for [M/H]=0 for 10 s yr (dotted line) and 5 x 10 9 yr (solid line); grainless models 
for [M/H] =-2 and t=10 Gyr (dashed line) The crosses, circles and squares correspond to 
observations from Dahn et al (1995), Monet et al (1992) and Dahn et al (1999), respectively. 
The indicated masses (in M ) correspond to the 5 Gyr [M/H]=0 isochrone. 

Fig. 13. — M^ vs (J-K) diagram for different ages and metallicites: [M/H]=0 for 10 8 yr 
(dotted lines) and 5 x 10 9 yr (solid lines); [M/H] =-2, t=10 Gyr (dashed line). The red 
curves on the right correspond to the DUSTY models for [M/H]=0. The blue curves on 
the left correspond to the COND models for [M/H]=0 (see text, §4.5.2). Filled circles and 
triangles on the isochrones indicate masses either in M or Mj (1 Mj « 1O~ 3 M ). Small 
green circles: MS stars from Leggett (1992) and Leggett et al (1996); black squares: Dahn 
et al (1999). Some identified BDs are also indicated (green diamonds). 

Fig. 14. — Deuterium-depletion curves as a function of age. The dashed lines correspond 
to different SSO masses, while the solid lines correspond to the 50% and 99% D-depletion 
limit, respectively. The inset displays the corresponding T e s(t). 
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Fig. 15. — Luminosity function in the V-band from the 5-pc sample (dash-line) (Kroupa, 
unpublished), the 8-pc sample (solid-line) (Reid et al 1999) and from the HST (dot-line) 
(Gould et al 1998) 



Fig. 16. — Mass function derived from the 5-pc LF (squares, last bin out of the figure) and 
the 8-pc LF (circles). Dash-line: dn/dm = 0.024m" 12 (see text) 
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